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B 1933+503, a dusty radio quasar at z > 2 : 
implications for blank field sub-mm surveys? 
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ABSTRACT 

We present a detailed mm-wave and optical study of the gravitational lens system 
B 1933+503, discovered by Sykes et al. (1998) in the radio. This object is probably 
the most complex lens system known, with 10 lensed components within a radius of one 
arcsecond. It is potentially important as a probe of the Hubble constant, although no 
optical counterpart has thus far been observed down to I — 24.2. We have obtained new 
sub-millimetre detections at 450 fim, 850 /im and 1350 ^m. We have also constrained 
the possible dust emission from the proposed foreground lensing galaxy using a if -band 
adaptive optics image and CO (5-4) measurements. A lensing model is constructed, 
taking the foreground elliptical galaxy at z = 0.755 as the lensing mass. From this we 
derive a scenario from which to model the sub- millimetre emission. Several arguments 
then point to the source in the B 1933+503 system lying above a redshift of 2. We 
speculate that unlensed relatives of this source may constitute a sizable fraction of the 
850 /im source counts. 

Key words: galaxies: active - galaxies: starburst - galaxies: individual: B 1933+503 
- gravitational lensing - cosmology: observations - infrared: galaxies 



1 INTRODUCTION 

The sub-millimeter waveband has recently become an in- 
valuable tool for investigation of the properties of high red- 
shift galaxies. An early result using the Sub-millimetre Com- 
mon User Bolometer Array (SCUBA) by Smail, Ivison and 
Blain (1997), showed that a much larger population of dusty, 
high redshift galaxies existed than previously thought. Since 
then, much work has gone into identifying just what types of 
galaxies make up this sub-mm bright population, and what 
mechanisms may be responsible for the rest-frame far-IR 
emission (Barger et al. 1998; Eales et al. 1998; Hughes et 
al. 1998; Ivison et al. 1998; Smail et al. 1998). An important 
application of SCUBA results is in constraining the global 
star formation rate, and a related unresolved issue is deter- 
mining the fraction of the sub-mm emission due to AGN 
activity, as opposed to a star-bursting process. 

The recently discovered radio object, B 1933+503, is 
among the most complex gravitational lens systems yet ob- 
served, with 10 lensed components within one square arc- 
second (Sykes et al. 1998). The existence of such a large 
number of images makes this a potentially important can- 
didate for determining the Hubble constant, through the 



measurement of relative time delays between image pairs. 
However Hubble Space Telescope V and I band images of 
B 1933+503 (also described in Sykes et al. 1998) reveal the 
lensing galaxy, but no apparent optical counterpart for the 
source down to about / = 24. Hence the source redshift has 
not yet been obtained. There are two scenarios which could 
account for the lack of optical emission: 1) the source object 
is intrinsically underluminous in the optical; or 2) the source 
or foreground lens is heavily dust enshrouded and thereby 
obscured at optical wavelengths. 

In order to characterise this object in the framework of 
dusty high redshift sources, and perhaps to shed some light 
on the source redshift, we obtained sub-millimetre contin- 
uum photometry at 450 /jm, 850 /im and 1350 fim using the 
SCUBA instrument on the James Clerk Maxwell Telescope. 
We have also indirectly constrained the dust obscuration in 
the lensing galaxy using an upper limit to CO line emission, 
and a high resolution, adaptive optics K band image. Us- 
ing the combined multi- wavelength data-set, we attempt to 
further constrain the properties of this system, and towards 
this end have also developed a new lensing model. 

The paper is organized as follows. We first describe the 
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Table 1. The observed properties of B 1933+503, and the foreground elliptical, from the radio to the optical waveband. None of the 
source quantities are corrected for the effects of gravitational lensing, which has a likely amplification factor of ~ 14 in the radio and 
perhaps 12 in the sub-mm and far-IR. 



Property 



Telescope 



B 1933+503 



Comment 



Distant AGN (source): 

a (J2000) 
S (J2000) 

Redshift estimate from lens model 

Flux density at: 

1.7 GHz 
5 GHz 
8.4 GHz 
15 GHz 

1350 /im 
850 /urn 
450 jim 

100 /im 
60 /tin 
25 (im 
12 (im 

K 

I 

V 

Foreground elliptical galaxy (lens): 

-Tap 
U ap 

12 CO (5-4) 
Redshift 



19 h 34 m 30?899 from VLA map component 4 (Sykes ct al. 1998) 
50° 25' 23'.' 22 

> 2.34 see Section 3.2 



VLA (D) 75.9 mjy All radio components (Sykes et al. 1998) 

VLA (D) 57.5 mjy All radio components (Sykes et al. 1998) 

VLA (D) 41.1 mjy All radio components (Sykes et al. 1998) 

VLA (D) 37.1 mjy All radio components (Sykes et al. 1998) 

JCMT 30±7mJy SCUBA photometry 

JCMT 24.1±2.6mJy SCUBA photometry 

JCMT 114±17mJy SCUBA photometry 

IRAS < 443 mjy XSCANPI limit 

IRAS < 136 mjy XSCANPI limit 

IRAS < 69 mjy XSCANPI limit 

IRAS < 69 mjy XSCANPI limit 

CFHT/AOB £ 20.3 Point source limit 

HST £ 24.2 Point source limit (Sykes ct al. 1998) 

HST £25.1 Point source limit 



CFHT/AOB 17.1 + 0.5 Aperture magnitude 

HST 20.6 ± 0.2 Aperture magnitude (Sykes ct al. 1998) 

HST > 22.5 Aperture magnitude limit (Sykes et al. 1998) 

JCMT <490mKkms" 1 Assuming 300 km s" 1 lincwidth 

Keck 0.755 ± 0.001 C. Fassnacht (private communication) 



new and previously existing observations of B 1933+503. We 
then describe a new lensing model using the foreground 
elliptical galaxy at z = 0.755 as the lensing mass. Various 
emission mechanisms are explored for the sub-mm flux, tak- 
ing into account the flat-spectrum radio components, dust, 
and the possible contamination of the signal by foreground 
emission in the z = 0.755 elliptical. From this we derive the 
most likely scenario for this object, and consider why it is 
not more readily observable in the optical and near-infrared 
bands. 



2 OBSERVATIONS 

B 1933+503 (hereafter 'B 1933') was discovered as part of 
the Cosmic Lens All-Sky Survey (Jackson et al. 1995), and 
has already been investigated in some detail at radio and 
optical wavelengths (see Sykes et al. 1998). We summarize 
the existing data, as well as our new observations, for both 
B 1933 and the proposed lensing galaxy in Table 1. We in- 
clude the radio data, as well as upper limits from IRAS, 
derived using the XSCANPI application provided by IPAC. 
These IRAS limits are the 90 per cent confidence limits for 
the flux per beam at the position of B 1933. 



2.1 Submillimeter observations 

Our observations were conducted with the SCUBA instru- 
ment (Holland et al. 1998) on the James Clerk Maxwell Tele- 
scope. On 1997 December 3, we operated the 91 element 
Short-wave array at 450 /im and the 37 element Long- wave 
array at 850 /im simultaneously in photometric mode, and 
also the single photometry pixel at 1350 /mi, giving half- 
power beam widths of 7.5, 14.7, and 21 arcsec respectively. 
At 450/850 pm five 'scans' of 900 seconds, consisting of 50 
integrations each were taken, for a total integration time 
of 1.25 hours. The central pixel of SCUBA was fixed on 
the source, B 1933. The 1350 /mi observation consisted of an 
additional integration time of one hour. In all cases, a 9- 
point jiggle pattern was employed to reduce the impact of 
pointing errors by averaging the source signal over a slightly 
larger area than the beam, resulting in greater photomet- 
ric accuracy. Whilst jiggling, the secondary was chopped at 
7.8125 Hz by 90 arcsec in azimuth, thus the other pixels ro- 
tate relative to the sky through the period of the integration, 
smearing out the effect of possible chopping on to other faint 
sources. Pointing was checked hourly on the blazar 2036-419 
and a sky-dip was performed between each 15 minute scan 
to measure the atmospheric opacity. The rms pointing er- 
rors were below 2 arcseconds, while the average atmospheric 
zenith opacities at 450 /mi 850 /mi and 1350 /mi were fairly 
stable with r being 0.51, 0.12 and 0.04 respectively. However, 
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b1 933+503 




Figure 1. CFHT ii"-band adaptive optics image of the 10 arc- 
sec X 10 arcsec field surrounding B 1933, smoothed with a 0.31 
arcsec boxcar filter (corresponding to the resolution). There is no 
indication of a central point source. 
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Figure 2. HST WFPC2 F814W (/-band) image, showing the 
same region as above, also smoothed to 0.31 arcsec resolution. 



there were some short time-scale variations, presumably due 
to water vapour pockets blowing over at high altitude, which 
caused some parts of the data-set to be noisier (see Borys, 
Chapman & Scott 1998 for more details). 

The data were reduced using both the Starlink pack- 
age SURF (Scuba User Reduction Facility, Jenness & Light- 
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Figure 3. Surface brightness profiles of the CFHT/AOB i\-band 
(top) and HST /-band (bottom) images, smoothed with a 0.31 
arcsec boxcar filter, and with an r 1 / 4 law fit to the central regions. 



foot 1998), and independently using our own routines (Borys 
et al. 1998). Spikes were first carefully rejected from the dou- 
ble difference data. The data were then corrected for atmo- 
spheric opacity and sky subtracted using the median of all 
the array pixels except for obviously bad pixels and the cen- 
tral pixels (the 1350 /im pixel currently has no provision for 
subtracting sky variations using the other wavelength pix- 
els). The data were then calibrated against Saturn and the 
compact HII region K3-50, which were also observed during 
the same observing shift. The two 850 /im and 1350 fxm cal- 
ibrations agreed with each other and also with the standard 
gains to within 10 per cent. However, at 450 /im, K3-50 is 
extended and variable, and is not a good calibration source, 
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while the Saturn 450 /im calibration agreed with standard 
gains to within 25 per cent. 

We also used the JCMT RxB3 345 GHz band hetero- 
dyne receiver on 1998 January 3, with the widest bandwidth 
possible (920 MHz) to observe the 12 CO(5-4) 575 GHz line, 
redshifted to 327.64 GHz. The intention was to check for 
a possible high gas mass from the lensing galaxy. No detec- 
tion was apparent in the data, which were binned to 38 MHz 
channels, with a resulting rms noise of around 4mK. We es- 
timate an upper limit of 0.49 K kms -1 for a Gaussian line 
with 300 km s -1 FWHM centred on the measured redshift 
z = 0.755 (see section 3.1). 

2.2 CFHT Adaptive Optics, If -band image 

We observed B 1933 with the Canada-France-Hawaii Tele- 
scope in March 1998 with the Adaptive Optics Bonnette 
(AOB, Rigaut et al. 1998) in place feeding the KIR near-IR 
camera. The detector has 1024 x 1024 pixels at a scale of 
0.34 arcsec per pixel to Nyquist sample the J-band under 
optimal observing conditions, for a field size of ~ 36 arc- 
sec. The average natural seeing throughout the observations 
was worse than median for the site at 0.86 arcsec, and the 
performance of the system was therefore not optimal, with 
fairly low Strehl ratios. The estimated resolution in the im- 
age near the V = 13.2 guide star is 0.22 arcsec correspond- 
ing to a physical scale of 890 h~ pc at z = 0.755 (assuming 
fio = 1 and with Ho = 100ft kms -1 Mpc -1 ). The adaptive 
optics corrected Point Spread Function (PSF) degrades in a 
predictable way with distance from the guiding source, as a 
function of guide star brightness, atmospheric quality, and 
wavelength of observation. Therefore the PSF for the po- 
sition on the detector of B 1933 was constructed from the 
guide star using a model describing the PSF variation over 
the observed field (Hutchings et al. 1998). B 1933 lies ~ 15 
arcsec from the guide star and the PSF has a FWHM of 0.31 
arcsec. 

Flux and PSF calibrations were performed using the 
UKIRT standard stars fsl3 and fs25 (Casali & Hawarden 
1992). Flat-field images were taken on the dome with the 
lamps turned on and off to account for the thermal glow of 
the telescope. The image was processed using an algorithm 
described in Chapman et al. (1998, in preparation), where 
the effects of a noisy, variable sky at 2.2 /im are reduced 
using a high signal to noise thermal image consisting of all 
the Jf-band sky data taken over a 5 night observing run. 
The signal to noise (~ 5) of the resulting image is too low 
to merit deconvolution techniques. The image is displayed, 
smoothed with a 0.31 arcsec boxcar filter, in Fig. 1, and 
the HST archive WFPC2 J-band image in Fig. 2. There is 
no obvious detection of the 4 lensed point sources of the 
radio core in either image, as discussed in the next section. 
Furthermore, even the foreground galaxy is not detected at 
V-band with HST. 

3 MODELLING 

3.1 The z = 0.755 elliptical 

Before considering the radio source itself, we examine the 
possibilities for the sub-millimeter emission to arise in the 



z = 0.755 foreground galaxy. The [I — K] colour of the galaxy 
is 3.5 ± 0.5, and [V — I] > 1.9, which is fairly normal 
for late-type galaxies at this magnitude and redshift range 
(e.g. Cowie et al. 1994). The colors are certainly consistent 
with some dust reddening though, with amplitude depend- 
ing on how red the V — I index is. We fit profiles to both 
the HST J-band and CFHT Jf-band images of this galaxy, 
as shown in Fig. 3. An r 1 / 4 profile is a good fit to the core 
in both cases, with no indication of a central cusp, i.e. no 
detection of the radio source. The HST J-band image shows 
unresolved knots to the northwest, extending along what ap- 
pears to be a spiral arm or tidal tail. The profile is distorted 
at this radius and is not well fit by a disc or bulge model. 
There is some evidence for elongation along PA = — 45° of 
the K-b&nd image, especially in the first 500 second integra- 
tion when the seeing was the best. This is consistent with the 
HST image, although the K isophotes are closer to circular. 

If we suppose that all the sub-mm emission comes from 
this galaxy instead of the radio source, we can calculate 
the dust mass required to produce this flux level (follow- 
ing e.g. Hughes, Dunlop & Rawlings 1997). The implied 
dust mass is M d ~ 2 x 10 8 /i -2 M Q , with an associated star 
formation rate of certainly >5OOM0yr~ 1 . This is a copi- 
ous amount of dust, exceeding that inferred for some ultra- 
luminous infrared galaxies (e.g. Eales and Edmunds 1996). 
For a fairly normal looking galaxy, such a large star forma- 
tion rate would be difficult to explain. In addition, we shall 
see in section 3.3 that the model fits to the spectral energy 
distribution (SED) are difficult to reconcile if the sub-mm 
emission is at 2 = 0.755. 

We can also use the non-detection of CO at the red- 
shift of the elliptical galaxy to estimate the maximum dust 
mass which could be present. Following Braine & Dupraz 
(1994), we can take the RMS level achieved at a particular 
binning scale, and estimate the la line intensity limit for 
lines of a given width. For 300 km s^ 1 FWHM lines we find 
that Jco <0.49K kms" 1 , or < MJykms" 1 for the JCMT 
at this frequency. The implied molecular hydrogen gas mass 
Mh 2 <>6x 10 9 ^ 2 M Q , taking CO (5-4) /CO (3-3) ~ 1 and us- 
ing a standard Mh 2 to Leo conversion (e.g. Barvainis et 
al. 1997). An accurate calculation of the mass of molecu- 
lar gas would require detailed modelling and multiple line 
constraints. However, if we take this value, together with a 
typical value for the gas to dust mass ratio, say ~ 500, then 
the dust mass is ^ 10 7 ft~ 2 Mo. Even taking the modelling 
uncertainties into account, this is still an order of magni- 
tude less than the 850 fim continuum measurement would 
imply. Uncertainties in the redshift and linewidth probably 
alos reduce this limit somewhat, but these effects still make 
it hard for the sub-mm flux to be coming from the elliptical 
galaxy. The redshift would have to be in error by around 
0.005, which is much larger than the estimated error, before 
our detection threshold would be much affected. Moreover, 
the dust mass depends on linewidth only as AJ d cx I 1 / 2 , 
which is a relatively weak dependence. Moreover, even as- 
suming a dust to gas ratio as small as ~ 100, the implied 
dust is still smaller than that implied by the 850 /im flux. 

Although the distorted, knotty morphology at J-band is 
suggestive of ongoing star formation, we can conclude from 
the above arguments that at most ~ 10 per cent of the 
sub-mm emission could be arising in the foreground ellip- 
tical. This does not imply that the foreground lens cannot 
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be a source of optical obscuration. Indeed a class of 'dusty 
lenses' has been discussed in the literature (Larkin et al. 
1994; Lawrence et al. 1995; Malhotra et al. 1996) where the 
background source is extremely reddened by small amounts 
of dust in the foreground lens. However, in at least one of 
these cases (MG1131+0456) the dusty lens hypothesis has 
been questioned after HST-NICMOS imaging of the source 
was obtained at //-band (Kochanek 1998). Without such a 
detection of our source, we are not able to rule out the pos- 
sibility of some reddening by the foreground lens. Therefore 
the obscuration of the optical source is currently an unknown 
combination of self-absorption and absorption through the 
lens. 



3.2 Gravitational lensing 

Nair (1998) modelled the lensing of B 1933 using a poten- 
tial derived from the shape of the HST /-band image of the 
foreground galaxy, and the positions of the detected radio 
sources. To accurately interpret the spectral energy distri- 
bution of B 1933, we first expanded on the Nair model using 
our new data, and then re-modelled the system to assess 
the most probable redshift for the source and the sub-mm 
to far-IR lensing amplification. We will give a more detailed 
analysis of the lens remodeling and the implications for time 
delays elsewhere (Lewis et al. in preparation). 

3.2.1 Physical parameters of the lens/source 

To estimate the source redshift, the crucial parameter to 
evaluate is the model mass-to-light ratio, T = r (M/ C) , 
which contains the quantity r — D s /Di- S , which is the ratio 
of angular size distances from the observer to the source and 
from the lens to the source. Here M and C are the mass and 
luminosity of the galaxy within some aperature. A constraint 
on r = T c /T , where T c is the expected mass-to-light ratio 
of a normal elliptical galaxy, can be used to constrain the 
redshift. 

We write T c = T p // evo i, where T p is the value appro- 
priate for present day ellipticals and / ovo i is the evolutionary 
correction accounting for the fact that galaxies were brighter 
in the past. We adopt a blue-light T p = (11.86 ± 0.5)h from 
van der Marel (1991) and note that the error is believed to 
be underestimated by about a factor of 3 due to the use 
of a fundamental plane relationship to determine distances. 
The evolutionary correction adopted is / evo i = 2.5 ± 0.5 
(— 1±0.2 magnitudes) from the models described in Pozzetti 
et al. (1996). 

The mass of the lensing galaxy is calculated out to a 
major-axis distance of 1 arcsec from the Nair model, and the 
luminosity is obtained from the reported photometry (Sykes 
et al. 1998) within a 2 arcsec x 1 arcsec aperature (which is 
slightly flatter than the model axial ratio of b/a — 0.59). 
The apparent blue magnitude of the lensing galaxy is calcu- 
lated from the expression in Lilly et al. (1996) which relates 
the observed /-magnitude to the blue magnitude using a K- 
correction derived from the CFRS data. The uncertainty in 
this estimate is dominated by the photometric errors (~ 20 
per cent). 

Two cosmological models, both with A = 0, were con- 
sidered: (i) a standard inflationary model with Qq = 1 and 



(ii) an open model with flo — 0.3. The results are indepen- 
dent of the Hubble constant and turn out not to be very sen- 
sitive to the cosmological model parameters. However, the 
derived source redshift is very sensitive to the parameters 
involved in the calculation of T c and T . Using the adopted 
mean values, we find r = 1.452 and 1.640 in cases (i) and 
(ii) respectively. These are close to the critical values cor- 
responding to an infinite source redshift and, formally, we 
find z s ~ 20 in either case. ^From a consideration of the 
uncertainties associated with the various parameters that 
enter into the calculation, we suggest that a realistic error is 
Ar/r ~ 50 per cent. Taking the upper limit on r, we obtain 
lower bounds of z s = 2.34 and 2.18 for the two cosmological 
models. A model SED is plotted in Fig. 5b for this redshift 
lower limit (for the Q.® = 1 model). 

There is a separate statistical argument about the 
source redshift, based simply on the lensing optical depth. 
Normally this is phrased in terms of the probability of find- 
ing the lens at some redshift, given a source redshift, but we 
can easily turn it around. Kochanek (1993) has discussed the 
gravitational lens statistics for isothermal spheres, showing, 
for example, that for a source at redshift z s , the most proba- 
ble location of the lens lies at half the proper motion distance 
to the source. For z\ — 0.755, and with a uniform prior dis- 
tribution on the redshift, we find a 95 per cent confidence 
lower limit of z B > 1.1, and a most likely redshift of 2.9 (for 
an Qo = 1 model) . While not particularly constraining, this 
lends some additional support to the above argument for the 
source being at z <; 2. 

3.2.2 Re-modelling of the lens system 

We employed a simple elliptical potential structure 
(Kochanek et al. 1989), and while the resulting mass profile 
is different from that used by Nair (1998), and we only used 
the observed positions as a constraint, the model recovers 
the image positions to similar accuracy. The total amplifi- 
cation of the core radio component in this model is 13.9, in 
reasonable agreement with the 14.63 found by Nair (1998). 

The degree to which the flux from a source is enhanced 
by the action of gravitational lensing is dependent on the 
scale-size of the source, leading to possibly pronounced dif- 
ferential magnification effects (e.g. Schneider & Weiss 1992). 
We model the sub-mm emission as coming from a dusty disc 
of ~ 200ft -1 parsccs in extent, lying perpendicular to the 
axis defined by the radio core and jet/lobes, and assuming 
f2o = 1 to convert lengths to angles. This corresponds to an 
extended obscuring torus possibly related to the small-scale 
(few pc) dusty torus of the standard unified AGN model 
(see for example the HST image of NGC4261 from Jaffe 
et al. 1998). Other sizes and geometries could of course be 
considered; we only present this as one possibility. 

Fig. 4 presents the results of our modelling, showing the 
source plane (left-hand panel) and the resulting image plane 
(right-hand panel). The darker shading represents both the 
steep and shallow spectrum radio emission (the central com- 
pact source is a little hard to see in the source plane). The 
field is 1.4 arcsec in extent and the image configuration can 
be compared directly to fig. 1 in Sykes et al. (1998). Lighter 
grey here represents the sub-mm emission from a dusty disc 
model as described above. The resulting far-infrared and 
sub-mm amplification (Air) is 11.9. 
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Figure 4. The result of our lens modelling for one particular example of a dusty disc. The left panel represents the source plane, while 
the right panel is the (observed) image plane. The solid lines denote caustics in the source plane and critical lines in the image plane. 
Each panel is 1.4 arcsec on a side, and is oriented to align with fig. 1 of Sykes et al. (1998). Darker shading represents both the shallow 
spectrum (compact source) and shallow spectrum (extended 'lobes') radio emission, while lighter grey shows a model of the far-infrared 
and sub-mm emission from a possible dusty torus. The resulting image configuration can be compared to the observed system. 



3.3 A dusty radio quasar 

All the available data for B 1933 from radio through optical 
are plotted as a spectral energy distribution in Fig. 5. The 
point source limiting magnitudes from the background noise 
in the K , I and V images are converted to flux units using 
the zero points described in Fukugita et al. (1995). Stickel 
et al. (1996) show that in many such galaxies, the near-IR 
can be substantially brighter than the optical. Our K-band 
image with a resolution of 0.31 arcsec FWHM allows us to 
put a 5o" limit of K — 21.4 on point source detections within 
the field. This is comparable to the limit of the HST image 
of I = 24.2 (Sykes et al. 1998) given a steep spectrum like 
the Stickel et al. (1996) quasar model. The 5<r point source 
sensitivity in the l/-band image is V = 25.1 

The sub-mm data points are used to constrain models 
consisting of a thick dusty torus, possibly associated with 
an ongoing burst of star formation, surrounding the cen- 
tral AGN (Granato et al. 1996). These models were shown 
to provide good fits to several high redshift ultraluminous 
systems, thereby throwing into contention the assumption 
that star formation is the dominant engine in objects such 
as IRAS F10214+4724 (see Granato et al. 1996, 1997). The 
standard grey-body emission from a dusty cloud is expected 
to dominate the sub-mm emission: F v oc v 13 B(y). We fix the 
dust emissivity and temperature at typical values for this 
type of object, f3 = 1.5, T = 80 K. A single parameter then 
describes the orientation of the dust torus with respect to 
the line of sight, and a nearly edge-on view (corresponding 
to 0°) is required to obscure the quasar continuum source 
enough to approach the upper limits of the K and / bands 
(see Fig. 5a, which is for the best estimate of the redshift, 



z = 2.9). None of the plotted curves actually falls under 
the near-IR limits. However, there is enough extra freedom 
in modelling the dust torus that one cannot conclude from 
this that extra obscuration is required. On the other hand 
small amounts of dust in the foreground lens can certainly 
be invoked to explain the lack of optical/near-IR counter- 
parts (Larkin et al. 1994; Lawrence et al. 1995; Malhotra et 
al. 1996). It is generally easier to fit the near-IR limits if the 
source is at higher redshift, as shown in Fig. 5b. However it 
must be stressed that HST near-IR observations are show- 
ing this hypothesis to be false in some objects (Kochanek 
1998). 

These model SEDs are difficult to fit to the sub- 
mm/FIR data if the redshift lies much below our esti- 
mated redshift limit (for Qo = 1) of 2 = 2.34. In Fig. 6 
we plot the sub-mm flux ratio 5 , 45o/S l 850 as a function of 
redshift for curves spanning a range of models of dusty 
star- forming/ AGN galaxies (similar to fig. 3 of Hughes et 
al. 1998), incorporating extreme parameters for optically 
thin cool dust emission (30 K < T d < 90 K; 1 < f3 < 2). We 
thereby place fairly robust upper and lower limits on the per- 
missible redshift of 1.3 < z < 3.5. However the actual geom- 
etry and emission characteristics of the source are unknown, 
and there are conceivable scenarios where these limits might 
not hold (e.g. if there were a much larger non-thermal con- 
tribution at 850 /im than assumed here). So again this ar- 
gument is suggestive only, although it is now the third such 
argument pointing to z > 2. 

The most statistically likely redshift is 2 = 2.9. How- 
ever, our lens model, with reasonable parameters for the 
foreground lens, is best fit with even higher values, and with 
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Figure 5. The spectrum of B 1933 from radio to optical, with 
an extrapolated radio power-law component, and a dusty AGN 
spectrum (see Granato ct al. 1996) fit to the sub-mm points: 
(a) 3 orientations of a dusty torus surrounding the central AGN 
assumed to be at z = 2.34 here - 0°, 45° and 90°; (b) the SED 
rcdshifted to 2 = 2.34 (limit from lens modelling), 2 = 2.9 (best 
guess redshift), and 2 = 4.5 (highest redshift consistent with this 
model SED). 
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Figure 6. 450 fim to 850 fimflux ratio for different models. The 
two solid lines span a range of reasonable emissivities and dust 
temperatures. 

arguments it seems fairly robust to assume z > 2 for the 
source. 

Although the 1350 fim point seems to imply that some 
radio power law component is still contributing to the emis- 
sion, the results for our dust emission model fit are not sub- 
stantially affected. If the non-thermal source contributes at 
850 /xm at all, a direct extrapolation from the radio shows 
it is minimal, and at most would affect the dust emissivity 
model by a factor of less than 0.2. 

We can estimate the dust mass using the opti- 
cally thin approximation as 1.0 x 10 8 A u ^Ii~ 2 Mq , with 
/3 = 1.5, T = 80 K, and Aib. the sub-mm lens amplification 
(and here using z s = 2.9 in an Qo = 1 universe). If we 
assume that some fraction /sf of the sub-mm emission is 
due to dust heated by young stars, rather than the cen- 
tral AGN continuum source, we can estimate a star forma- 
tion rate of 400 x /sf h~ 2 . Calzetti, Kinney & Storchi- 
Bergmann (1994) have shown that the amount of obscura- 
tion is highly dependent on the distribution of the dust. The 
fact that the optical/near-IR counterparts to the bright ra- 
dio nucleus are not observed can easily be explained by a 
combination of dust obscuration in B 1933 itself, and dust 
in the foreground lens. 



a lower limit to the redshift of z = 2.34. It is also clear from 
this dust torus model that redshifts in excess of z = 4.5 can 
remain consistent with the sub-mm and far-IR data points, 
especially if the dust temperature is increased (see Fig. 5b). 
There is a reasonable lower limit on the redshift from our 
sub-mm data, and certainly by 2 = 0.755, the redshift of the 
foreground lens, they are not even consistent with the IRAS 
upper limits. From the weight of these individual redshift 



4 DISCUSSION 

Given that the radio properties of B 1933 suggest that it is 
a radio-loud elliptical, it is reasonable to assume that the 
sub-mm emission results from a roughly edge-on dusty disc 
surrounding the core of the galaxy and heated by the AGN. 
However, if the AGN phase of massive galaxies has any con- 
nection with enhanced star formation rate, then it is cer- 
tainly possible that some fraction of the sub-mm emission 
results from dust heated by young massive stars. The radio 
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luminosity is, however, not a good measure of the AGN sub- 
mm flux, and so we cannot ascertain what fraction of the 
dust is heated by star formation rather than AGN activity 
without a spectrum of the source. 

An HST-NICMOS program (H.-W. Rix, private com- 
munication) will image B 1933 at //-band in the near future, 
and may reveal the source, including some morphological in- 
formation. This should add a strong constraint to the lensing 
predictions. Our work predicts that the //-band magnitude 
must be greater than 23, so that it may be within the 5a 
detection threshold of the HST observation for a reasonable 
exposure time. 

Although the beamsizes of the sub-mm measurements 
are fairly large, substantially larger than all the lensed com- 
ponents, the average SCUBA source density has been shown 
to be ~ 2000 degree" 2 (e.g. Smail et al. 1997, Hughes et 
al. 1998, Eales et al. 1998) at the depth of - 2.5 mjy RMS 
that we reach at 850 fj,m (the most sensitive of our observa- 
tional windows) , and it is unlikely that another dusty source 
contributes to the sub-mm flux. We expect only one 24 mjy 
source per 0.1 degree 2 and therefore the chance for a spuri- 
ous detection of an unrelated source is negligible. 

Blain has (1996) modelled the effects of lensing on the 
distribution of sub-mm sources. A large fraction of the galax- 
ies detected at 850 fim flux densities in the range 0.1-10 mjy 
could be gravitationally lensed, implying that lensing in the 
sub-mm may be a very important consideration for high 
redshift source counts. Systems such as B 1933 at ~2 mjy 
intrinsic flux density, as well as the Cloverleaf quasar and 
IRASFW2 14+4724 (Barvainis et al. 1995) are known lensed 
sub-mm sources, which are relatively extreme objects in the 
amount that they are amplified. We have then caught a 
rather rare view of a somewhat average dusty radio galaxy at 
an early epoch, thanks to an excellent chance alignment with 
a foreground galaxy. This high redshift object appears to be 
very much a scaled down version of its more extreme cousins 
(such as SMM 02399-0136 (Ivison et al. 1998), 8C1435+635 
(Ivison et al. 1998), or APM 08279+5255 (Lewis et al. 1998), 
which are otherwise the only objects bright enough to study 
in detail. 

It is interesting to note that if this object were unlensed, 
or only lensed by the potential of a rich cluster such as the 
sources discovered in cluster surveys (Smail et al. 1997, 1998; 
Chapman et al. in preparation), then the sub-mm emission 
at 850 and 450 fim would still be detectable by SCUBA in a 
deep blank field integration, but the 1350 /im emission would 
be undetectable with SCUBA, and the radio components 
would be hard to detect in a VLA integration, such as that 
used for the discovery of B1933 with a la RMS of ~0.4mJy 
(Sykes et al. 1998). The optical identification would then be 
very difficult, given that the true source might have / > 27. 
There would be no way to assess the AGN contribution to 
the source, and the object would likely be identified as a star 
forming galaxy with a fairly large SFR of > lOOMoyr" 1 . 
We hypothesize that a substantial fraction of the faint dusty 
galaxies discovered in sub-mm surveys could harbour AGN 
as their main power source, affecting the inferred comoving 
SFR at high z. Although optical surveys have shown that 
high redshift quasar systems are relatively rare, this could 
of course be affected by dust obscuration, as in the case 
of B 1933. Indeed, serendipitous discovery of heavily dust- 
enshrouded young QSOs at z ~ 2 (e.g. H 167, Q 0059-2735, 



Egami et al. 1996) implies a possibly sizable population of 
optically obscured AGN at high-z. 

Another class of dusty objects, the EROs - extremely 
red objects with [R - K] > 6 (e.g. Hu & Ridgway 1994), 
have also recently been detected in the sub-mm (Cimatti et 
al. 1998). It is still unclear whether there is any AGN con- 
tribution in these objects; B 1933 might be another example 
of this class of object. 

The contribution of AGN to counts of SCUBA-bright 
galaxies is difficult to estimate. This is particularly true since 
the volume surveyed is essentially unknown (because the 
redshift range is currently unknown) . In addition there is the 
possible extra effect of lensing amplification. However, the 
case of B 1933 certainly leads us to be cautious in interpret- 
ting all SCUBA detections as highly star-forming galaxies. 



5 CONCLUSIONS 

We have presented new sub-mm and near-IR observations 
of the complex gravitational lens system, B 1933+503. Our 
estimate of the total dust mass is likely to be correct to 
within the errors in our lensing magnification estimate, since 
the slope of the SED shows that optically thin dust is the 
most viable explanation for the emission in this wavelength 
region. The fact that the optical/near-IR counterparts to the 
bright radio nucleus are not observed can easily be explained 
via dust obscuration in B 1933 itself, or by invoking a 'dusty 
lens' (as in Lawrence et al. 1995). The large sub-mm flux 
shows that this object is certainly not underluminous, as 
originally proposed. 

The results show the power of sub-millimeter observa- 
tions for detection and characterization of obscured galaxies 
in the distant universe. Investigation of this highly amplified 
AGN leads us to wonder what fraction of sources detected 
in 'blank' SCUBA fields could be AGN, or lensed, or even 
both. 
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